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<f> meson production is studied by the NA49 Collaboration in central Pb+Pb collisions at 20A, 
30A, 40A, 80A, and 158A GeV beam energy. The data are compared with measurements at lower 
and higher energies and with microscopic and thermal models. The energy dependence of yields 
and spectral distributions is compatible with the assumption that partonic degrees of freedom set 
in at low SPS energies. 



I. INTRODUCTION 



The production of strange particles is considered one of the key observables for understanding the reaction mecha- 
nisms in ultrarelativistic heavy-ion collisions. Enhanced strangeness production with respect to proton-proton colli- 
sions was originally proposed as a signature of the transition to a deconfined state of quarks and gluons during the initial 
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stage of the reactions The enhancement was predicted to arise from gluon fragmentation into quark-antiquark 
pairs which is believed to have a significantly lower threshold than strange-antistrangc hadron pair production chan- 
nels. Indeed, it has been observed [2, H] that the ratio of the number of produced kaons to that of pions is higher by 
a factor of about 2 in central S + S and Pb + Pb reactions than that in p + p collisions at the top energy available 
at the CERN Super Proton Synchroton (SPS). 

Statistical hadron gas models have been successfully employed to describe the measured particle yields at various 
collision energies [3, H, 0, Q ■ The fact that the hadronic final state of the collision resembles a hadron gas in 
chemical equilibrium has been interpreted as a consequence of the hadronization process [9|] or as a result of a fast 
hadronic equilibration process involving multiparticle collisions (Io| . In this hadron gas picture, enhanced production 
of strange particles in collisions of large nuclei arises as a consequence of the increased reaction volume, relaxing 
the influence of strangeness conservation [ill ]. Technically, this requires the application of the canonical ensemble 
to small collision systems, while for larger volumes such as those encountered in central collisions of heavy ions, the 
grand-canonical approximation is valid. It has been shown that this "canonical strangeness suppression" also applies 
to a partonic system fl2l ]. 

In addition to this volume effect, the strange particle phase space appears to be undersaturated in elementary 
interactions. The deviation of the str ange p article yields from a hadron gas in full equilibrium was parametrized by a 
strangeness undersaturation factor 7s [3 [l3[ . The additional suppression becomes much weaker in heavy-ion collisions. 
However, fits to the hadron multiplicities in full phase space are still unsatisfactory when not taking into account 75 
Q. A possible interpretation is that the total amount of strangeness available for hadronization is determined in a 
prehadronic stage of the collision. A change in 7s between p + p and A + A would then reflect the difference in the 
initial conditions of the respective fireballs. 

The hadron gas model was extended to describe the energy dependence of produced hadron multiplicities by a 
smooth paramctrization of the fit parameters T and [Ib-, determined at energies available at the BNL Alternating 
Gradient Synchroton (AGS), SPS, and BNL Relativistic Heavy Ion Collider (RHIC), as a function of collision en- 
ergy [l4j . However, this extended model failed to reproduce the detailed features of the energy dependence of relative 
strangeness production measured by NA49 in its energy scan program. In particular, the sharp maximum at around 
30A GeV beam energy [ll| [l6| could not be described. The same holds true for microscopic reaction models such as 
UrQMD [l7|. On the other hand, this feature was predicted as a consequence of the onset of a phase transition to a 
deconfined state at the respective beam energy [l8| . 

In this context, it is certainly interesting to investigate specific strangeness-carrying hadrons. Among these, the 
(j) meson is of particular interest because of its ss valence quark composition. In a purely hadronic scenario, being 
strangeness-neutral, it should not be sensitive to hadrochemical effects related to strangeness. If on the other hand, 
the amount of available strange quarks is determined in a partonic stage of the collision, the <f> is expected to react 
more sensitively than singly strange particles. In particular, one would expect the <fi meson yield to be suppressed 
by 7^ with respect to equilibrium. Analogously, the canonical suppression mechanism in small systems should have a 
stronger effect on the </>, leading to a larger relative enhancement in Pb + Pb collisions with respect to p + p reactions 
than observed for kaons. 

In the evolution of the fireball after hadronization, <f> mesons can be both formed by kaon coalescence and destroyed 
by rescattering. In addition, when decaying inside the fireball, the daughter particles can rcscatter, leading to a loss 
of signal in the invariant mass peak of the respective decay channel. This is more likely to happen for slow <p mesons, 
which spend more time in the fireball. Thus the effect could lead to a depletion of the 4> meson yield at low pt in 
central nucleus- nucleus collisions [l9| . 

Theoretical investigations have suggested that the properties of the <p meson might be modified in a dense hadronic 
medium. In particular, a decrease of its mass of the order of 10 MeV 1 [2(| and an increase of its width by a factor 
of 2-3 [21 1 were predicted. So far, there is only one experimental claim for a broadening of the width in p + Cu 
collisions 22| . 

In an earlier publication (23|, we reported on (f> production at top SPS energy, where we found the <fi enhanced by 
a factor of about 3, compared to minimum bias p + p collisions at the same beam energy. Meanwhile, the <p meson 



was measured at the same energy by the NA50 [24|, NA45 [25j, and NA60 [26| experiments. At the AGS, data on 
4> production were obtained by the E917 Collaboration in Au + Au collisions at Pbcam = 11-7 A GeV (^/snn = 4.88 
GeV) in a restricted rapidity range [j/ c . m . — 0.4, y c m ] [27j . At the RHIC, the STAR Collaboration measured the cj> 
meson at ^/snn = 130 and ^Jsnn = 200 GeV at midrapidity [H, |2j|. For the latter energy, data are also available 
from the PHENIX experiment [3fl|. 

In this article, we report on (f> production in central Pb+Pb collisions at five different beam energies from 20A 



1 For better readability, we use natural units, i.e., c = 1, throughout this article. 
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TABLE I: Characteristics of the data sets employed in the analysis. The mean numbers of wounded nucleons (N w ) were 
obtained by Glauber model calculations. 
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(N w ) 


N events 


Momentum range 


(A GeV) 


(GeV) 












(GeV) 


20 


6.3 


1.88 


2002 


7.2% 


349 ± 1 ± 5 


352 309 


2.0-23.0 


30 


7.6 


2.08 


2002 


7.2% 


349 ± 1 ± 5 


368 662 


2.0-27.0 


40 


8.8 


2.22 


1999 


7.2% 


349 ± 1 ± 5 


586 768 


2.0-27.0 


80 


12.3 


2.57 


2000 


7.2% 


349 ± 1 ± 5 


300 992 


2.0-32.0 


158 


17.3 


2.91 


1996 


5.0% 


362 ± 1 ± 5 


345 543 


3.5-35.0 



to 158A GeV. Together with the data obtained at the AGS and the RHIC, our findings enable the study of energy 
dependence of </> production over a large range of collision energies. 

II. EXPERIMENT 

The NA49 experiment at CERN is based on a fixed-target hadron spectrometer using heavy-ion beams from the SPS 
accelerator. Rs main components are four large-volume time projection chambers for charged-particle tracking, two 
of which operate inside the magnetic field of two superconducting magnets, thus providing an excellent momentum 
measurement. Two larger main time projection chambers (MTPCs) are placed downstream, outside of the field, 
and enable particle identification by the measurement of the specific energy loss in the detector gas. The particle 
identification capabilities are enhanced by a timc-of-flight (TOF) scintillator system behind the MTPCs, albeit in a 
restricted geometrical acceptance. 

A thin lead foil with 1% interaction probability for Pb nuclei was used as a target. For the different runs, the 
magnetic field was scaled proportionally to the beam energies in order to have similar acceptance in the cm. system. 

The centrality of the reactions was determined from the energy deposited by the beam spectators in the zero-degree 
calorimeter, placed 20 m downstream of the target. By setting an upper limit on this energy, the online central trigger 
selected the 7.2% most central collisions at 20A-80A GeV and the 10% most central collisions at 158A GeV. The 
latter data set was restricted to 5% centrality in the offline analysis. The corresponding mean numbers of wounded 
nucleons were obtained by Glauber model calculations (see Table H|). Details of the experimental apparatus can be 
found in Ref. [p]. 

III. DATA ANALYSIS 
A. Event and track selection 

Offline quality criteria were applied to the events selected by the online centrality trigger to suppress nontarget 
interactions, pileup, and incorrectly reconstructed events. The cut variables include the position and x 2 of the 
reconstructed vertex and the track multiplicity. For the central data sets used in this analysis, however, the impact 
of these quality cuts is marginal; only about 1% of all events were rejected. Table U shows the event statistics used in 
the analysis for the five data sets. 

The analysis was restricted to tracks reconstructed in the MTPCs which could be assigned to the primary vertex. 
A minimal track length of 2 m out of the maximal 4 m in the MTPCs was required to suppress ghost or split tracks 
and to ensure a good resolution in dE/dx. Detailed studies including reconstruction of simulated tracks embedded 
into real raw data events showed that for such a selection of tracks, losses due to track reconstruction and high track 
density are negligible. 

B. Selection of kaon candidates 

NA49 observes the <j> meson through its hadronic decay into charged kaons. To reduce the large contribution of 
pions and protons to the combinatorial background, kaon candidates were selected based on their specific energy loss 
dE/dx in the MTPCs. The mean dE/dx of pions, kaons, and (anti-)protons was determined from TOF-identificd 
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FIG. 1: (Color online) dE/dx parametrization for the data set at 80^4 GeV. (a) Mean dE/dx as function of /3y determined for 
TOF-identified pions, kaons, and protons; (b) dE/dx resolution as function of momentum, obtained from the deconvolution of 
the energy loss spectra into the contributions of ir + , K + , and p. 



particles in the acceptance of the time-of- flight detectors and parametrized as a function of f3j as shown in Fig.QJa). 
This allowed one to extend the momentum range for the identification from the TOF acceptance to higher momenta. 
The lower momentum limit was given by either the MTPC acceptance or the crossing of the Bethe-Bloch curves of 
pions and kaons. The momentum limits for the different data sets are summarized in Table HJ 

Fixing the mean dE/dx of kaons and protons to this parametrization. the resolution was obtained by unfolding the 
energy-loss spectra in momentum bins into the Gaussian contributions of the particle species (p, K, ix, and e). The 
resolution is about 4% and has a slight momentum dependence which was again parametrized [Fig.QJb)]. 

Kaon candidates were selected by a momentum-dependent dE/dx window around the expectation value, the size of 
which was chosen to optimize the <j) signal quality. In addition, the window had to be symmetric and large enough to 
minimize the sensitivity to the errors in the determination of the dE/dx expectation value and resolution. A window 
of ±1.5cr was found to be the best choice. This selection contains 87% of all kaons, giving an efficiency of 75% for 
the pair. The fraction of true kaons within the selected candidate track sample varies between 40% and 60%. 



C. Extraction of raw yields 



The 4> signal was obtained by calculating the invariant mass of all combinations of positive and negative kaon 
candidates in one event. To reconstruct the combinatorial background of uncorrelated pairs, candidates from different 
events were combined. The mixed-event spectrum was subtracted from the same-event spectrum after normalization 
to the same number of pairs [32j ■ Figure [5] shows the background-subtracted invariant-mass spectra in the total 
forward acceptance for different collision energies. In all cases, clear signals are observed at the expected position. 

While the subtracted spectrum is flat on the right side of the signal, a depletion is observed between the peak and 
the threshold. As a possible source of this undershoot, the correlation of kaons stemming from different <p mesons 
has been discussed in Ref. [32| . In our case, it was shown by simulation that this effect is small thanks to the large 
acceptance of the NA49 MTPCs. Another possible source of the distortion is the reflection of other resonances, 
e.g., A — * ir~p, into the K + K~ spectrum by misidentification of pions and protons, as discussed in detail in 
Ref. [HI. This effect was shown to be present in our previous analysis of another data set [23-], where the dE/dx 
resolution was significantly worse. However, all such resonances would distort the spectrum over a broad range above 
threshold, which can be excluded by the observed flatness at higher masses. This conclusion is further strengthened 
by the observation that the depletion does not vanish when applying a stricter dE/dx cut on the kaons. Hence, the 
undershoot is likely to originate from a true correlation of kaon pairs. Simulations show that it can be explained 
by final state strong interaction of kaons (34|. This is demonstrated in Fig. [3] by showing the K + K~ correlation 
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FIG. 2: (Color online) K + K~ invariant-mass spectra after subtraction of the combinatorial background in the forward rapidity 
hemisphere for the five different beam momenta. The full lines show the Breit-Wigner fits to the signals as described in the 
text. The bin size is 4 MeV for 20A and 30^4 GeV and 2 MeV for the other beam energies. 



function in gi nv = y/ (pi — P2) 2 — (E\ — E2) 2 [Fig.[3fa)] and in mi„ v [Fig.^b)]. While the repulsive interaction causes 
a depletion in m mv , the stronger attractive Coulomb effect is squeezed into 0.8 MeV above threshold and is thus 
hardly seen. In combination with the steeply rising unsubtracted mi nv distribution, this depletion can easily account 
for the deficit observed in the subtracted spectrum. 

To correct this effect quantitatively by simulation is difficult and would moreover be model dependent. As the 
narrow signal is easily distinguished from the broad residual background, we accounted for the depletion by fitting 
a straight line in the vici nity of the peak. For the description of the signal itself, we used a relativistic p-wave 
Breit-Wigner distribution [35[ of the form 

dN mT(m) 

dm (m 2 — m 2 ) 2 + m 2 ] T 2 {m) 

with the mass-dependent width 

T{ m ) = 2vJ^] (2) 
\<W q 2 + % 

where q := <J jTO 2 — m 2 K and qo := <J \rriQ — m 2 K . This distribution was folded with a Gaussian representing the 

invariant-mass resolution a m of the spectrometer. Since in general, mass resolution and width cannot be determined 
separately, we fixed the width to its book value To = 4.26 MeV (36[, leaving mo, cr m , a normalization and two 
parameters for the linear background as free parameters for the fit, which was performed in the mass range 994- 
1050 MeV. It was checked by simulations that this procedure gives the correct values for position, width and integral 
of the distribution. As Fig. [^demonstrates, the fit gives a good description of the signal. The numerical values of the 
fitted parameters are listed in Table HT1 

To obtain longitudinal and transverse spectra, the signal was extracted in rapidity and in p t bins, respectively, in 
the same way as in the total acceptance. Generally, the limited statistics prevented a simultaneous division into y-pt 
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FIG. 3: K + K correlation function close to threshold in (a) qi nv and (b) mi nv [34 



bins. Thus, transverse momentum spectra could only be derived averaged over rapidity. To reduce the number of free 
fit parameters, to and a m were fixed for the fits in the phase space bins to the values obtained from the signal in 
the total acceptance. For the 158A GeV data set, where the statistics in the signal allowed to do so, we checked that 
leaving these parameters free did not significantly alter the results. In particular, no significant dependence of too or 
<7 m on rapidity or p t was observed. 

Since the straight-line background is only an approximation for the residual background in the vicinity of the 
signal, the stability of the fit against the variation of the fit region was checked. The parameters too and a m show no 
significant dependence. The variation of the normalization constants, which determine the fit integral, is in all bins 
far below the statistical error returned by the fit procedure. We conclude that the latter properly takes into account 
the possible variations of the baseline. 

The raw yields in the phase-space bins were obtained by integrating the fit function from threshold up to too+30Fo ~ 
1.148 MeV. This mass cutoff is somehow arbitrary; the corresponding integral varies by about 3% for cutoff values 
from too + 10Fo to infinity. We take this as a systematic uncertainty due to the mass cutoff. Using alternatively a 
(analytically intcgrablc) nonrclativistic Lorentz distribution for the fit does not change the integral by more than 1%. 



D. Geometrical acceptance 

The geometrical acceptance of the NA49 detector for the decay <fi — > K + K~ was obtained double-differentially in y 
and p t (integrated over azimuth) by GEANT simulations of the (f> decay including in-flight decay of the kaon daughters, 
assuming an azimuthally flat emission and isotropic decay. The resulting acceptance is shown in Fig. [4] for 20 A and 
158 A GeV. While the upper momentum limit for the daughter candidates restricts the acceptance at forward rapidity 
for the top SPS energy, at lower beam energies there is lack of acceptance near midrapidity because of the lower 
momentum limit for the daughter tracks and the increased losses due to in-flight decay for low-momentum kaons. 

As the acceptance is a function of y and p t , the proper correction factor for a given extended phase-space bin 
(integrated either over y or p t ) as used in the analysis is the mean acceptance 

J dydp t a(y,p t ) f(y,p t ) 
— s 

Jdydpt f{y,p t ) 
s 

where S denotes the region in the y,pr plane, a(y,pt) the acceptance probability averaged over the azimuthal angle, 
and fijjiPt) the differential <p meson yield. For the rapidity distributions, the differential yields have in addition to be 
extrapolated to the full p t range. The extrapolation factor, however, is small (<5%) due to the large p t range covered 
by the experiment. 
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FIG. 4: Geometrical acceptance probability for <j> — > K + K including kaon decay in flight for (a) 20^4 GeV and (b) 158A GeV. 



Both the acceptance correction and the extrapolation to full pt require the knowledge of the y and pt dependence 
of 4> meson yields, which leads to an iterative procedure (see Sec. IIII E|) . 



E. Spectra and yields 



Apart from the differential acceptance correction, the raw yields obtained from the fit to the invariant-mass spectra 
were corrected for the branching ratio </> — > K + K~ (49.1%) and the efficiency of kaon dE/dx selection (75% for the 
pair), and normalized to the number of collisions. These global correction factors are common for all bins in phase 
space and for all beam energies. 

The transverse spectra are fitted by the thermal ansatz 

^oc Pt e—/ T , (4) 
dp t 



where the transverse mass mt = \Jfn^ -Vp\- The distributions in rapidity were parametrized by a single Gaussian 



dn 
dy 



oc e 



(5) 



As the parameters T and a y must be obtained by the analysis itself, an iterative procedure was employed. Starting 
from some reasonable parameter values, the acceptance correction was calculated according to Eq. ([3]), assuming 
factorization of the emission function f(y,Pt) into the transverse and longitudinal distributions (U) and i.e., 
independence of T on rapidity. The corrected yields in the p t and y bins were then fitted with the distributions ^ 
and ([5]), respectively, obtaining new values for T and a y which serve as input for the next iteration. Convergence of 
the method was reached after three to five steps. It was checked that the final results do not depend on the choice of 
start values for the parameters. 

After the final step of the iteration, the yields in full phase space were obtained by summing up the measured 
yields in the rapidity distributions and numerically extrapolating Eq. ([5]) to the full rapidity range. In a similar way, 
the quantities (pt), (mt), and a y were determined. The midrapidity yield dn/dy was obtained directly from the fit 
function. 

As demonstrated later in Fig. [3 the Gaussian parametrization gives a satisfactory description of the rapidity 
distribution for all data sets. However, because of the lack of midrapidity data points at the lower beam energies, an 
ambiguity for the extrapolation to full phase space arises. To check the sensitivity of the results to the assumed shape 
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of the rapidity distribution, we alternatively parametrized the latter by the sum of two Gaussian functions displaced 
symmetrically around midrapidity by a shift a: 



dn 

— oc e y 
ay 



(0) 



The width of this distribution will be characterized by its rms value. Total yield (</>) , midrapidity yield dn$ j dy and 
rms y were calculated for both parametrizations ([5]) and ©. The final values listed in Tables fVl and I VII were calculated 
as the mean of the results of the two methods; their differences enter the systematic errors. 



F. Statistical and systematic errors 



Statistical errors in the raw differential <f> meson yields originate from the statistical bin-by-bin errors in the same- 
event and mixed-event invariant mass spectra, which were found to be in good approximation Poissonian and un- 
corrected between mass bins. Then, the statistical errors in the event-mix subtracted invariant-mass spectrum was 
calculated as HH 

of = n ,i + k 2 n em>i , (7) 

where no,j is the number of entries in mass bin i in the same-event spectrum, n em j the same number in the mixed- 
event spectrum, and k the normalization constant for the event mix. These errors were propagated toward the raw 
differential yields by the least-squares fit of the Breit-Wigner distribution to the signal peak. 

The acceptance calculation was performed with sufficiently high statistics such that the relative statistical error 
of the differential acceptance is below 1% and thus far below the uncertainty in the raw yields over the entire y, pt 
region used for the analysis. Finally, the errors in the acceptance-corrected differential yields are propagated through 
the least-square fits to the spectra to obtain the statistical uncertainties in the spectral parameters and the integrated 
quantities. 

Systematic uncertainties in the uncorrected yields arise from the approximation of the residual background in the 
invariant-mass spectra as a straight line. This approximation is only valid in a limited mass range around the signal 
peak. Thus, the stability of the results of the Breit-Wigner fit against the variation of the fit range was checked. 
We found no significant dependence of the parameters mo and a m ; the variation of the normalization constant, 
determining the fit integral, was in all y and pt bins found to be smaller than the statistical error. 

Another source of systematic error arises from the dE/dx selection of kaon candidates. Uncertainties in the 
paramctrization of the mean kaon dE/dx and the resolution result in systematic deviations of the efficiency cor- 
rection from its true value. To estimate this error, the analysis was repeated for different widths of the dE/dx 
selection window around the kaon expectation value, applying the respective efficiency correction. This error was 
found to be the dominating one; for most raw yields, it is comparable to or slightly larger than the statistical one. 

Imperfect detector description in the simulation leads to systematic uncertainties in the acceptance correction. To 
reduce possible errors, the analysis was restricted to phase-space regions where the acceptance is above 1%. The 
remaining error was estimated by repeating the analysis with varying acceptance conditions (minimal track length 
in the MTPCs). It was in all cases found to be much smaller than the error originating from the kaon selection by 
dE/dx. 

As the spectral parameters enter the acceptance correction through Eq. ([3]) , their uncertainties add to the systematic 
errors of the corrected yields. This was accounted for by determining the range of acceptance values allowed by the 
errors in T and a y . In addition, for the rapidity bins close to beam rapidity, a possible deviation of the slope parameter 
by 50 MeV from its averaged value was taken into account in the acceptance correction. The resulting error, however, 
is small thanks to the large and approximately uniform p t acceptance. 

The systematic errors in the corrected differential yields were assumed to be independent and added in quadrature. 
They were propagated to the respective errors in the spectral parameters by repeating the fit of Eqs. ((U)-© with 
statistical and systematic errors added and comparing the resulting errors to those obtained from the fit with statistical 
errors only. 

For the determination of the averaged quantities (</>), (p r ), (mt}, and rmSj,, the summation of the measured differ- 
ential yields as well as extrapolation to full phase space are required. The systematic errors of these observables were 
determined from the errors of the differential yields and the uncertainties in the spectral shapes. 
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TABLE II: Approximate number of detected <f> mesons S, background-to-signal ratio B/S, signal-to-noise ratio SNR, position 
of the signal peak mo, and invariant-mass resolution <r m . The latter two were obtained by a Breit-Wigner fit to the signal peak 
(see text). The width was fixed to its literature value 4.26 MeV. S and B were calculated in a window of ± 4 MeV around the 
peak. The quoted errors are statistical only. 



.Pbcam 


S 


B/S 


SNR 


m 


dm 


(A GeV) 








(MeV) 


(MeV) 


20 


6 500 


70 


9.4 


1018.8 ±0.6 


2.6 ±0.9 


30 


16 500 


104 


12.5 


1018.4 ±0.5 


2.5 ± 1.3 


40 


37 000 


53 


26.2 


1018.9 ±0.2 


2.1 ±0.3 


SO 


55 000 


30 


42.5 


1019.1 ±0.1 


1.1 ±0.1 


158 


180 000 


72 


49.4 


1019.0 ±0.1 


1.8 ±0.1 



IV. RESULTS 
A. Line shape 

Tabic |TT] summarizes the parameters obtained from the invariant-mass signals in the total acceptance. The signal 
quality decreases when going to lower beam energy because of both the reduced <j> meson yield and the reduced 
acceptance due to the increased in-flight decay probability for the daughter kaons. At all five energies, the fitted 
peak position is slightly below the literature value of 1019.43 MeV [36[. We investigated the effect of an error in the 
normalization of the magnetic field used for momentum determination in the reconstruction chain and found that a 
bias of 1% in the magnetic field is needed to explain the observed shift. This is slightly above the momentum scale 
uncertainty deduced from a precision study of the K® signal. We thus cannot exclude that the deviation of the peak 
position is due to experimental effects. 

The widths of the mass peaks obtained from the fits are consistent with those obtained from a full detector simulation 
and reconstruction. Their slight increase toward lower beam energies can be understood as the increasing influence of 
multiple scattering on lower momentum tracks. For the signal at 158A GeV, we fitted simultaneously width and mass 
resolution and obtained Tq = (4.41 ± 0.61) MeV, a m = (1.81 ± 0.26) MeV, i.e., no deviation from the free-particle 
width. Thus, within experimental uncertainties, we do not observe indications for a mass shift or a broadening of the 
(j) meson. 

The observation that the mass and width of the <j> meson agree with the Particle Data Group values is in line with 
the results of AGS and RHIC experiments [U [H, 0, [3(| • It should be noted that because of the long lifetime of the 
</> meson (r = 46 fm), only a fraction decays inside the fireball. Thus, only a part of the </> mesons can be expected to 
be influenced by the surrounding medium. 

B. Transverse momentum spectra 

The transverse momentum spectra obtained for the five beam energies are shown in Fig. [SJ numerical data are given 
in Table Ifffl In all cases, the thermal distribution Q gives a good description of the data; the fit parameters are 
summarized in Table lIVl At top SPS energy with the best signal quality, a modest deviation from the fit function is 
indicated by the % 2 /ndf of 1.5. A slight curvature of the transverse mass spectrum at this energy, as expected from a 
hydrodynamical expansion scenario, is visible for this energy in Fig. [5Ja). For the other energies, no deviations from 
pure exponential behavior can be seen within the experimental uncertainties. 

The transverse momentum spectrum can be also characterized by its first moment or the average transverse mass. 
These parameters were calculated from the measured data points and extrapolated to full p t using the exponential fit 
function. As the extrapolation contributes only marginally because of the large pt coverage, (pt) and (m t ) — mo are 
largely independent of the spectral shape. Their values are also listed in Tabic [TVl 

The assumption of the slope parameter being independent of y could be checked for 158A GeV, where statistics 
allowed us to extract transverse spectra in four different rapidity bins. The resulting slope parameters are shown in 
Fig. [B](b) ■ Within the measured rapidity range, we observe no significant change of the slope parameter with y. Using 
the y-dependent slope parameters for correcting the rapidity distribution had no sizable effect on the results. 

The spectrum obtained for 158 A GeV agrees with that from an earlier publication [23| of the NA49 experiment, 
which was based on the analysis of an older data set at the same beam energy. For comparison, the previously 
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FIG. 5: (Color online) <f> transverse momentum spectra integrated over the rapidity intervals given in Table ITVl The full lines 
show the fits of thermal distributions Q. The squared symbols denote previously published results (2^. Only statistical errors 
are shown. 
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FIG. 6: (Color online) (a) cA transverse mass spectra integrated over the rapidity intervals given in Table HVl The exponential 
fits indicated by the full lines correspond to the fits shown in Fig. [S] The spectra for different beam energies are scaled for 
better visibility. Only statistical errors are shown. The data at 158^4 GeV are compared with previously published results of 
NA49 [23j | and CERES [25|. (b) Slope parameter as function of rapidity at 158^4 GeV. The values agree within errors with that 
obtained from the {/-integrated pt spectrum, the latter indicated with its standard deviation by the shaded bar. 
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TABLE III: Differential (f> meson yields in the p t (left) and y (right) distributions. Data in the pt bins are integrated over the 
rapidity ranges given in Table ITVl The errors are statistical. 



Pt (GeV) 


dn/ (dydpt) (GeV~ 


- 1 \ 

1 


y 


dn/dy 






E beam = 20A GeV 






0.0-0.4 


0.382 ± 0.074 




0.2-0.6 


1.043 ±0.250 


0.4-0.8 


0.528 ± 0.097 




0.6-1.0 


0.536 ± 0.077 


0.8-1.2 


0.257 ±0.054 




1.0-1.4 


0.159 ±0.033 


1.2-1.6 


0.079 ± 0.030 




1.4-1.8 


0.032 ± 0.017 


1.6-2.0 


0.033 ±0.015 












E beam = 30A GeV 






0.0-0.3 


0.231 ±0.051 




0.3-0.6 


0.735 ±0.194 


0.3-0.6 


0.578 ± 0.079 




0.6-0.9 


0.651 ± 0.090 


0.9-1.2 


0.386 ± 0.079 




0.9-1.2 


0.456 ± 0.052 


1.2-1.5 


0.257 ± 0.050 




1.2-1.5 


0.193 ± 0.036 


1.5-1.8 


0.070 ±0.019 




1.5-1.8 


0.097 ± 0.029 






Ebeam = 40,4 GeV 






0.0-0.2 


0.185 ±0.035 




0.3-0.6 


1.067 ±0.108 


0.2-0.4 


0.668 ± 0.052 




0.6-0.9 


0.756 ± 0.059 


0.4-0.6 


0.780 ± 0.064 




0.9-1.2 


0.611 ±0.038 


0.6-0.8 


0.625 ± 0.075 




1.2-1.5 


0.348 ± 0.028 


0.8-1.0 


0.569 ± 0.073 




1.5-1.8 


0.188 ±0.023 


1.0-1.2 


0.413 ±0.059 








1.2-1.4 


0.275 ± 0.040 








1.4-1.6 


0.081 ±0.028 








1.6-1.8 


0.086 ±0.019 








1.8-2.0 


0.057 ±0.014 












E he am = 80A GeV 






0.0-0.2 


0.337 ±0.031 




-0.3-0.0 


1.591 ±0.304 


0.2-0.4 


0.886 ±0.051 




0.0-0.3 


1.474 ±0.138 


0.4-0.6 


1.148 ±0.057 




0.3-0.6 


1.258 ± 0.086 


0.6-0.8 


0.996 ± 0.056 




0.6-0.9 


1.351 ±0.062 


0.8-1.0 


0.861 ±0.052 




0.9-1.2 


1.041 ± 0.049 


1.0-1.2 


0.517 ±0.048 




1.2-1.5 


0.718 ± 0.043 


1.2-1.4 


0.344 ± 0.045 




1.5-1.8 


0.408 ± 0.037 


1.4-1.6 


0.173 ±0.040 




1.8-2.1 


0.197 ±0.040 






Ebeam = 158A GeV 






0.0-0.2 


0.582 ± 0.053 




0.0-0.2 


2.557 ±0.166 


0.2-0.4 


1.275 ±0.086 




0.2-0.4 


2.386 ±0.121 


0.4-0.6 


1.924 ± 0.098 




0.4-0.6 


2.229 ± 0.098 


0.6-0.8 


2.016 ±0.099 




0.6-0.8 


2.202 ± 0.089 


0.8-1.0 


1.778 ± 0.092 




0.8-1.0 


1.974 ± 0.090 


1.0-1.2 


1.339 ± 0.080 




1.0-1.2 


1.816 ±0.094 


1.2-1.4 


0.956 ± 0.067 




1.2-1.4 


1.636 ±0.105 


1.4-1.6 


0.567 ±0.055 




1.4-1.6 


1.528 ±0.126 


1.6-1.8 


0.370 ± 0.044 




1.6-1.8 


1.125 ±0.171 


1.8-2.0 


0.200 ± 0.034 









published data are shown by the square symbols in Figs. [5je) and[6ja). There is agreement with the results of the 
CERES experiment in both decay channels 4> — ► K + K~ and (j) — > e + e~ [251 ]. as also demonstrated in Fig. EJa). The 
data disagree with the spectrum measured by the NA50 experiment in the di-muon decay channel — > fi + /i~ , where 
a significantly smaller slope was obtained [24[ . 
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TABLE IV: Rapidity range (in cm. system), pt range, slope parameter T, % 2 P er degree of freedom, average p t , and average 
mt for the transverse momentum spectra. T and \ 2 are results from the fit of Eq. ((4]) to the spectrum; (pt) and (rru) — mo were 
obtained by summation over the data points and extrapolation to full p t using the fit function. The first error is statistical, 
the second one systematic. 



Pbcam (A GeV) 


y range 


Pt range (GeV) 


T (MeV) 


X >df 


(pt) (MeV) 


(m t ) - m [MeV] 


20 


0.0-1.8 


0.0-2.0 


196.8 ± 19.5 ±20.2 


1.06/3 


650.9 ± 34.2 ± 40.2 


229.5 ±20.1 ±23.6 


30 


0.0-1.8 


0.0-1.8 


237.4 ± 17.8 ±22.9 


2.03/4 


738.9 ± 28.3 ± 46.3 


284.6 ± 17.3 ±28.4 


40 


0.0-1.5 


0.0-2.0 


244.6 ±9.0 ±5.8 


12.42/8 


763.4 ± 15.8 ± 14.3 


297.8 ±10.0 ±9.2 


80 


0.0-1.7 


0.0-1.6 


239.8 ±8.3 ±10.9 


3.48/6 


756.4 ±11.5 ±22.5 


292.6 ±7.6 ±15.3 


158 


0.0-1.0 


0.0-2.0 


298.7 ± 6.6 ± 10.6 


12.06/8 


883.5 ±9.9 ±21.3 


378.3 ±6.7 ±15.2 



TABLE V: Parameters of the single-Gauss fit ((5} and the double-Gauss fit ([6| to the rapidity distributions. The RMS was 
calculated from the data points and extrapolated to the full rapidity range using the average of the two parametrizations. The 
first error is statistical, the second one systematic. 



Pbeam 




Xi/ndf 


a-. 


a 


Xi/ndf 


rmSj 


(A GeV) 














20 


0.572 ±0.037 ±0.030 


0.042/2 


0.425 ± 0.026 ± 0.022 


0.425 


0.75/2 


0.582 ± 0.031 ± 0.040 


30 


0.752 ±0.047 ±0.057 


2.02/3 


0.538 ± 0.028 ± 0.032 


0.538 


1.03/3 


0.769 ± 0.030 ± 0.062 


40 


0.863 ± 0.033 ± 0.042 


3.14/4 


0.696 ±0.118 ±0.036 


0.487 ±0.149 ±0.051 


3.10/3 


0.852 ± 0.015 ± 0.038 


80 


1.016 ± 0.028 ±0.033 


17.55/6 


0.658 ± 0.035 ± 0.043 


0.682 ± 0.029 ± 0.043 


4.12/5 


0.974 ± 0.024 ± 0.074 


158 


1.451 ± 0.086 ±0.012 


2.36/7 








1.444 ± 0.021 ± 0.054 



C. Rapidity distributions and yields 

Figure [7] shows the rapidity distributions, which for all five energies arc in good agreement with both the single- 
Gaussian and the double-Gaussian paramctrization (see curves). Numerical data are given in Table ED For the data 
sets at 2(L4 and 3CL4 GeV, due to the low number of data points, the double-Gaussian fit was constrained to a = <j y 
as suggested by the data at 40^4 and 80A GeV. 

Only at 80 A GeV is the complete forward hemisphere covered. At 158^4 GeV, large rapidities are not measured 
because of the upper momentum cut on the secondary kaons. Since kaons below 2 GeV laboratory momentum 
cannot be reliably identified by dE/dx because of the crossing of the Bethe-Bloch curves, no signal could be extracted 
at midrapidity for the lower three beam energies. The uncertainties in the extrapolation toward midrapidity is 
demonstrated by the difference of the two parametrizations. It adds to the systematic error of the total yield and, in 
particular, to that of dn/dy at midrapidity. Table fVl lists the parameters obtained by the two fit functions, respectively 

Alternatively, the rapidity distributions can be characterized by their second moments in a model-independent 
fashion. The root mean square of the distributions was calculated from the measured data and extrapolated to the 
full rapidity range using the parametrizations (O and ((6|) . The average of the two results is listed in Table [Vj 

Total yields were obtained by summation of the data points in the rapidity spectra and extrapolation to the full 
rapidity range by the average of the fit functions. The midrapidity yield dn/dy was calculated analytically from the 
average of the fit functions. For the determination of statistical and systematic errors, the correlation of the spectral 
parameters were properly taken into account. The results for the mean multiplicity of mesons (<$) and for the 
midrapidity yield dn^jdy are listed in Table PVTl 

All results obtained at 158^4 GeV are consistent within statistical errors with NA49 results published earlier [23| . 
which were obtained from a data set taken in 1995 (squared symbols in Fig.[7£e). The main difference of the two data 
sets is an improved dE/dx resolution, resulting in a reduced pion contamination of the kaon candidate sample. The 
cleaner kaon identification reduces the distortions in the background-subtracted invariant-mass spectrum induced by 
resonances with a pion as decay daughter (33|, thus leading to a smaller systematic error of the Brcit-Wigner fit to 
the spectrum. We thus prefer to use the newly obtained results at 158A GeV for the discussion. 
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FIG. 7: (Color online) (j> rapidity distributions. The solid points refer to measured data, the open points are reflected at 
midrapidity. The full lines show the parametrization by a single Gaussian ([5}, the dashed lines that by the sum of two 
Gaussians ©. The data at 158^4 GeV are compared with previously published results of NA49 [13] and CERES [2f|. Only 
statistical errors are shown. 



TABLE VI: Total (f> multiplicity (<j>) and midrapidity yield dn^/dy calculated from the rapidity distributions of Fig. [7] The 
first error is statistical, the second one systematic. 

Pbcam (A GeV) (#) dn /rfy(y c . m .) 

20 1.89 ±0.31 ±0.22 1.17 ± 0.23 ± 0.38 

30 1.84 ±0.22 ±0.29 0.94 ± 0.13 ± 0.30 

40 2.55 ±0.17 ±0.19 1.16 ± 0.16 ± 0.14 

80 4.04 ±0.19 ±0.31 1.52 ± 0.11 ± 0.22 

158 8.46 ±0.38 ±0.33 2.44 ± 0.10 ± 0.08 



V. DISCUSSION 



The enhancement of relative strangeness production in heavy-ion collisions with respect to proton-proton reactions 
is a well-known fact. In an earlier publication (23j . the enhancement factor for the (j> meson at top SPS energy was 
found to be 3.0 ± 0.7, thus larger than for kaons and A, but smaller than for multistrange hyperons. We calculate the 
<f> enhancement by normalizing the measured (f> meson yield in A ± A by the number of wounded nucleon pairs and 
dividing by the corresponding yield in p + p. For the lower beam energies, no reference measurements in elementary 
collisions are available. Here, we employ a parametrisation of the <f> excitation function in p + p collisions as described 
in Ref. (27[. For top SPS energy and RHIC, the </> meson yield measured in p + p [23|, [29| was used. Figure [5] shows 
the resulting enhancement factor 

2 {4>)a+a . . 

N w (<p)p+ P 

as a function of energy per nucleon pair. The measurement of the E917 Collaboration at AGS (Pbeam = 11-7A GeV) 
was extrapolated to full phase space assuming the same rapidity distribution as for K~ as suggested by the authors [27j • 
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FIG. 8: (Color online) tj> enhancement factor E$ [see Eq. (JSj) ] as function of energy per nucleon pair. Data from the AGS [2?| 
and the SPS refer to multiplicities in full phase space, data from the RHIC [2^, [3(| to midrapidity yields. The shaded boxes 
represent the systematic errors. 



In the AGS/SPS energy region, the value of E$ lies between 3 and 4, and within our experimental uncertainties we 
find no systematic variation here. At RHIC energies, the enhancement appears to be lower, significantly so, should 
the PHENIX result be validated. It should be noted, however, that the RHIC values were derived from midrapidity 
data while at lower energies phase-space integrated yields were used. 

In the context of statistical models, the enhancement of strangeness production can be interpreted as a result of the 
release of suppression due to strangeness conservation when going from small (p + p) to large (central A + A) systems. 
Technically, this is reflected in the application of the canonical ensemble for small systems, while large systems can 
be described by the grand-canonical ensemble. In this picture, a smaller enhancement at RHIC energies points to 
the fact that at such high energies, strangeness is produced with sufficient abundance for the canonical suppression 
to be relaxed even in p + p collisions. However, in a purely hadronic picture, canonical suppression does not act 
on the <j) meson because it is a strangeness-neutral hadron. Enhanced <p production can thus be attributed either 
to enhanced strangeness production in a partonic stage of the collision or to the coalescence of kaons which suffer 
canonical suppression also in a hadronic scenario. 

The hadrochemical models have been extended not only to fit hadron multiplicities for a given reaction but also to 
describe the energy dependence of particle yield ratios by a smooth variation of the relevant parameters T and \ib 
with collision energy 

mm- Here ' thc 

energy dependence of temperature and baryochemical potential is obtained 
by a parametrization of the values for T and fis obtained from fits to particle yield ratios at various collision energies. 
The model reproduces many yield ratios of the bulk hadrons; however, this does not hold for the 4> meson, as shown 
in Fig. Ht^a), where the measured excitation function of the {(/))/ (w) ratio [{tt) = 1.5((7r + ) + (tt~))} is compared with 
the model prediction. The relative cf> meson yields at thc SPS are overpredictcd by factors of up to 2. The situation 
remains essentially unchanged when midrapidity ratios are considered instead of integrated yields [Fig.[5Jb)]. At the 
RHIC, there is a large experimental ambiguity as a result of thc different results on <f> production obtained by the 
STAR and PHENIX experiments HHH. 

A better description of the data is obtained if a strangeness saturation parameter j s is allowed. The corresponding 
model predictions [8[ for the <j> multiplicity, resulting from a fit to the hadron abundances at 11.7A, 30A, 40A, 80A, 
and 158 A GeV, are compared with the data in Fig. [10] (solid points). Note that this model does not provide a 
continuous description of thc energy dependence; thc points arc only connected to guide the eye. Thc agreement with 
the measurements at the higher SPS energies is very good. Thc successful application of the saturation parameter 
7 S on the strangeness-neutral <f> meson for Pbeam > 40Al GeV again suggests that the strangeness content at chemical 
freeze-out is determined on a partonic level for these energies. 

Final state interactions after chemical freeze-out could change the equilibrium <j) yield and spectra. In particular, 
scattering of the daughter kaons with other produced hadrons would lead to a loss of the </> signal in the experimentally 
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FIG. 9: (Color online) {4>)/{tt) ratio (a) in full phase space and (b) at midrapidity as function of energy per nucleon pair [{n) — 
1.5((7r + ) + (7r~))]. The CERES data point [25| was displaced horizontally for visibility. Note that the CERES measurement is at 
y ~ 2/c.m. — 0.5. The full line shows the predictions of the extended hadron gas model (HGM) with strangeness equilibration [37f ] . 
the dashed curves those obtained with UrQMD 1.3 [13]. The shaded boxes represent the systematic errors. 




FIG. 10: (Color online) <j> multiplicity in central A + A collisions as function of energy per nucleon pair. The solid points denote 
the results of the statistical hadronization model (SHM) which allows a deviation from strangeness equilibrium They are 
connected by the solid line to guide the eye. The dotted curve shows the cf> yield predicted by the UrQMD 1.3 model [13]. The 
shaded boxes represent the systematic errors. 



observed decay channel, predominantly at small rapidities and low values of pt- Such a loss is not expected in the 
leptonic decay modes, since electrons or muons will leave the fireball without interaction. A comparison of the 
measured mt spectrum via the K + K~ and e + e~ decay channels [see Fig. EJa)] indicates that the effect cannot be 
large. To study the effect on the total yield, we used the string- hadronic transport model UrQMD [l7j . It was found 
that only about 8% of the decayed cf> mesons are lost for detection due to rescattering of their daughter particles, 
independent of collision energy. Similar results have been obtained with the RQMD model [l9(. The effect is thus 
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FIG. 11: (Color online) (a) Widths of the rapidity distributions of tt~ , K + , K~ , and <f> in central Pb + Pb collisions at SPS 
energies as function of beam rapidity [3Sl . l39fl . The dashed lines are to guide the eye. The open star denotes the (j> rapidity 
width measured in p + p collisions . (b) Widths of the <f> rapidity distributions in central Pb+Pb collisions compared with 
the expectations in a kaon coalescence picture [Eq. JjJ)]. The shaded boxes represent the systematic errors (shown only for <j> 
mesons) . 



not sufficient to account for the deviation of the relative 4> multiplicities from their equilibrium values. 

On the other hand, (j> mesons can be produced by KK scattering. In fact, kaon coalescence is the dominant (« 
70%) production mechanism for the </> in UrQMD, again for all investigated collision systems. As shown by the dotted 
curve in Fig. 1 1 01 the model gives a reasonable description of the cf> meson yields at lower energies, whereas it starts to 
deviate from the measurements at intermediate SPS energies. The discrepancy with data is more pronounced when 
studying the {(j>) / (n) ratio (Fig. [5]) because UrQMD overestimates the pion yields at SPS energies by about 30%. 

The hypothesis that the 4> meson is produced predominantly by kaon coalescence can be tested by comparing the 
4> and kaon distributions in phase space. Figure [TTl a) shows the width of the <f> rapidity distribution as a function of 
beam rapidity at SPS energies, together with that measured for 7r~, K + , and K~ [38l. l39j| . The (f> meson width does 
not fit into the systematics observed for the other particle species but increases much faster with energy. While at 
20 A GeV, the <f> rapidity distribution is narrower than that of K~ , we find it at top SPS energy comparable to the 
pions. In addition, at 158 A GeV it is much larger in central Pb + Pb collisions than measured in p + p collisions at 
the same energy J23j|, a feature which is not observed for other particle species. 

In the kaon coalescence picture, there would be a tendency for the <f> rapidity distribution to be narrower than those 
of the kaons. In an ideal case, neglecting correlations, 

— — _L 

2 2 ' 2 ' \ ) 

a l a K+ K~ 

where the distributions were approximated by Gaussians. As shown in Fig. lllf b). the <f> data rule out kaon coalescence 
as dominant formation mechanism for beam energies above 30A GeV. Only at 2QA GeV, the observed rapidity widths 
are consistent with the coalescence picture. As mentioned before, this would also explain the <p enhancement at low 
energies, where a transient deconfined state is not expected. 

The observation that models based on a purely hadronic reaction scenario have serious problems in describing 
relative strangeness production in the upper SPS energy range is not unique to the </> meson but holds for kaons and 
other strange particles, too. It has been related to the onset of deconfincment at around 30A GeV as predicted by 



the statistical model of the early stage [181 ]. A striking experimental evidence is the narrow maximum in the K + /n + 
ratio at this energy [3Sl. [39l] . A similar structure is, within experimental errors, not observed for the <fi meson (Fig. [9]); 
instead, the energy dependence of the relative 4> meson yield resembles that of the K~ . This can be understood since 
the K + yield is in good approximation proportional to the total strangeness production, which is not the case for K~ 
and 4> because a large, energy-dependent fraction of s quarks is carried by hypcrons. 



16 



> 

CD 



450 



600 



400 



350 - 



300 



250 - 



200 - 



150 - 



100 



A E9 1 7 

■ NA49 

NA50 

* CERES (K + K~) 

• STAE 

O PHENIX 



10 



(a) 



2 

10 

Vs NN [GeV] 



> 

CD 

2 550 r 
o 

f= 500 - 

I 

A 

^"450 - 

V : 

400 r 

350 r 

300 r 

250 r 

200 r 

150 - 



A E917 

■ NA49 

NA50 

• STAR 

O PHENI) 



ti ■ 



10 



Hi 



(b) 



10 

Vs NN [GeV] 



FIG. 12: (Color online) (a) Inverse slope parameter T and (b) average transverse mass (m t ) — mo of the (f> meson in central 
A + A collisions as function of energy per nucleon pair. The data from E917 [27| were averaged over the measured rapidity 
interval (see Table [[Vj- Results from NA50 [H| and RHIC (HlHHj were obtained at midrapidity, the result from CERES 
at y = —0.71. Data from NA49 are integrated over rapidity. The PHENIX data point was slightly displaced horizontally for 
visibility. For the NA49 data, {mt) was calculated from the transverse momentum spectra using an exponential extrapolation 
to full p t . For the other data sets, it was derived analytically from the exponential fit function. The shaded boxes represent 
the systematic errors. 



The energy dependences of both the inverse slope parameter and the mean transverse mass of the <j) meson arc shown 
in Fig. [TH The transverse mass spectra of the <j> are well described by exponential fits [see Fig. HJfa)]; consequently, 
the two parameters show a similar behavior. Over the energy range AGS-SPS-RHIC, there is an overall tendency 
for both parameters to increase. However, a constancy of the values in the lower SPS energy range, as has been 
observed for pions, kaons, and protons (39j — a fact interpreted as being consistent with a mixed partonic/hadronic 
phase [4(| — cannot be excluded. 



VI. SUMMARY 



We have presented new data on <fi production in central Pb+Pb collisions obtained by the NA49 experiment at 
20A, 30A, A0A, 80A, and 158A GeV beam energies. No indications of medium modifications of the <f> meson mass 
or width were observed. The energy dependence of the production characteristics was studied by comparing them 
with measurements at AGS and RHIC energies. We find that at low SPS energy, the data can be understood in a 
hadronic reaction scenario; while at higher energies, hadronic models fail to reproduce the data. A statistical hadron 
gas model with undersaturation of strangeness gives a good description of the measured yields. This suggests that (j> 
production is ruled by partonic degrees of freedom, consistent with the previously found indications for the onset of 
deconfinement at lower SPS energy. 
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